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ABSTRACT 

We present Advanced Camera for Surveys observations of MACS Jl 149.5+2223, an X-ray luminous galaxy 
cluster at z=0.544 discovered by the Massive Cluster Survey. The data reveal at least seven multiply-imaged 
galaxies, three of which we have confirmed spectroscopically. One of these is a spectacular face-on spiral 
galaxy at z = 1 .491, the four images of which are gravitationally magnified by 8 <J p, <^ 23. We identify this as 

anL* {Mb — -20.7), disk-dominated (B/T <^ 0.5) galaxy, forming stars at ~ 6 M Q yr _1 . We use a robust sample 
of multiply-imaged galaxies to constrain a parameterized model of the cluster mass distribution. In addition 
to the main cluster dark matter halo and the bright cluster galaxies, our best model includes three galaxy- 
group-sized halos. The relative probability of this model is P(A r haio = 4)/P(A r haio < 4) > 10 12 where AWo is the 
number of cluster/group-scale halos. In terms of sheer number of merging cluster/group-scale components, 
this is the most complex strong-lensing cluster core studied to date. The total cluster mass and fraction of 
that mass associated with substructures within R < 500 kpc, are measured to be M tot = (6.7 ± 0.4) x 10 14 M Q 
and / su b= 0.25 ±0.12 respectively. Our model also rules out recent claims of a flat density profile at 7<r 
confidence, thus highlighting the critical importance of spectroscopic redshifts of multiply-imaged galaxies 
when modeling strong lensing clusters. Overall our results attest to the efficiency of X-ray selection in finding 
the most powerful cluster lenses, including complicated merging systems. 

Subject headings: cosmology: observations — galaxies: clusters: individual (MACS J 1 149.5+2223) — galax- 
ies: evolution — gravitational lensing 



1. INTRODUCTION 

Strong gravitational lensing by galaxy clusters is a well- 
established probe of the extragalactic Universe, offering a 
magnified view of high redshift galaxies that would oth- 
erwise be beyond th e reach of present day telescopes 
(e.g. iFranx etail 119971: lEllis et all l200lt iKneib et all l200l 
ISmail et alJl2007tlSwinbank et alJl2007t iBradley et al.ll2008l) . 
The detailed cluster mass models required to intepret such ob- 
servations contain a wealth of information about the mass and 
structure of cluster cores , against which theoretical predic- 
tions can be tested (e. g. [Smith et afj 120051: IComerford et al.l 
2006: lSandetli]|2008l) . 

The majority of spectroscopically confir med strong lensin, 
clusters are at low redshi ft, i.e. z^O-3 (e.g. lKneib et al.lll99' 
Broadhurst 1it~aT1 12005b ISmith etal.1 l2005t iLimousin etaD 
20071: iRichard et al.l2009l) . In contrast, only four spectroscop- 



ically confirmed strong lensing clusters are known at z>0.5 
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dGladders et alJ 12001 IfaadaetaLl ppl iBorvs et ail 12001 
ISharon et alT]2()05t ISwinbank et all 120071: lOfek et al] 120081) . 
The Massive Cluster Survey (MACS: lEbeling et al.ll2001l) of- 
fers an unprecedented opportunity to expand the available 
sample of strong lensing clusters at 0.3 < z <^ 0.7. We present 
new results from this search: MACS J 1 149.5+2223 (hereafter 
MACS Jl 149; 11:49:34.3 +22:23:42.5 [J2000]) at z= 0.544, 
one of a complete su bsample of 12 MACS clusters at z > 0.5 
(Ebeli ng et al]|2007l) . 

In ^2] we describe the data; modeling and results are pre- 
sented in §|3] and summarized in ^4] We assume Hq= 
70kms- 1 Mpc" 1 , n M = 0.3 and n A = 0.7; at z = 0.544 1" corre- 
sponds to 6.35 kpc. All uncertainties and upper/lower limits 
are stated and/or plotted at 95% confidence. 

2. OBSERVATIONS AND DATA ANALYSIS 

MACS J 1149 was observed on 2004, April 22 with the 
Advanced Camera for Surveys (ACS) on-board the Hubble 
Space Telescope (HST) 10 for 4.5 ksec and 4.6ksec through the 
F555W and F814W filters respectively (GO: 9722, PI: Ebel- 
ing). These data were reduced using standard MULTIDRIZ- 
ZLE routines onto a 0.03"/ pixel grid. The reduced data re- 
veal a striking multiply-imaged disk galaxy comprising three 
tangential images (Al. 1/2/3; Fig. [TJ, and an additional image 
(A1.4) likely caused by part of the galaxy's disk lying adja- 
cent to the radial caustic in the source plane. A blue image 
pair (A2.1/2) also lies ~ 30" South East of the BCG, with its 
counter-image (A2.3) ~ 15" to the South West. A further five 
triply-imaged galaxies are identified based on their distinctive 

10 Based in part on observations with the NASA/ESA Hubble Space Tele- 
scope obtained at the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. 
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FIG. 1 . — TOP LEFT - V555//8i4-band color picture showing the full ACS field o f view . Red contours show the luminosity density of cluster galaxies - cluster 
halos were centered on the 5 luminous structures labeled A-E in the lens model ( j|3.2t . The white boxes marked #1, #2 and #3 show the regions displayed in 
more detail in the three numbered panels. PANEL #1 - The central ~ 80" X 80" of the cluster showing the multiple image systems discussed in the text. The 
cyan (outer), magenta and yellow curves show the z = 1.491, z = 1.894 and z = 2.497 tangential critical curves respectively. The inner cyan curve shows the radial 
critical curve for z = 1.491. PANEL #2 - A faint triply-imaged galaxy next to a cluster galaxy within a group of galaxies ~ 50" North of the BCG. PANEL #3 - 
Two candidate triply-imaged systems adjacent to a bright cluster galaxy ~ 100" North of the BCG; A6.4 marks the location of a possible fourth image of A6. 
BOTTOM ROW - Zoom into the four images of Al; morphol ogic al feat ures used to constrain the lens model are marked by numbered white circles. The black 
crosses and circle in the A1.2/A1.1 panels are discussed in § j|3.2l & [3~5l N is up E is left in all panels. 



colors and morphologies: A3. 1/2/3 and A4. 1/2/3 both lie be- 
tween the BCG and a dense group of cluster ellipticals ~ 50" 
to the NNW; A5 is embedded in the halo of a cluster ellipti- 
cal in the same group to the NNW of the BCG; A6. 1/2/3 and 
A7. 1/2/3 surround a bright elliptical galaxy ~ 100" North of 
the BCG. Numerous other faint blue background galaxies can 
be seen through the cluster core, however the lack of concor- 



dant colors and morphologies preclude a reliable identifica- 
tion of them as being multiply-imaged at this time. 

MACS J 1 149 was obser ved with the Lo w Resolution Imag- 
ing Spectrograph (LRIS: lQke et al.ll 19951) on the Keck-1 10-m 
telescope 11 on 2004, March 28 and 2005, March 6 employ- 

The W. M. Keck Observatory is operated as a scientific partnership 
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FIG. 2. — Stacked one-dimensional and two-dimensional spectra of Al.l, 
A1.2 and A1.3. The 1" wide slit on each target sampled both the central 
bulge and the disk of each image. 



ing a single multi slit mask per run. In 2004, we used the 
400/3400 grism, D560 dichroic, and 400/8500 grating cen- 
tered at 7200 A. A total integration of 10.8 ksec yielded the 
redshift of A2.1/2 as z= 1.894 via detection of Lyman-a in 
emission at 3519 A (A2.3 was also confirmed at z = 1.894 
in 2005). This redshift was used to constrain a preliminary 
lens model from which the redshift of Al was predicted to 
be z ~ 1.5 ±0.1. In March 2005 we centered the 831/8200 
grating at 8750 A to search for redshifted [oil] from Al. 
The resulting 1-hr spectrum resolved the [oil] doublet at 
A bs = 9282.5 A in all three of Al. 1/2/3, placing this galaxy at 
Z = 1 .4906 ±0.0002 (Fig.E). The 2004 observations also iden- 
tified A3.1 and A3. 2 at z = 2.497 via detection of the Lyman 
break, plus interstellar Silicon and Carbon absorption lines. 

3. MODELING AND RESULTS 

3.1. Luminosity Density Map 

To gain an initial view of the structure of the clus- 
ter mass distribution we selected galaxies with spectro- 
scopic or photometric redshifts be tween 0.5 <. z ^ 0.6 from 
Ma et al.'s (in prep.; see also iMa et alj 120081) catalog. 
This catalog is based on spectroscopic observations with 
DEIMOS on Keck II (yielding 217 cluster members), plus 
panoramic B/V /R c /I c /z' -band imaging with Suprime-CAM 
on the Subaru 8.2-m telescope and similar M*-band data from 
MegaPrime on the CFHT 3.6-m telescope. The resulting lu- 
minosity density map i s adaptively smooth ed to 3er signifi- 
cance using ASMOOTH dEbeling et al.ll2.Q06h and is presented 
in Fig. Q] It reveals five luminous structures (labeled A-E), 
suggesting that the underlying distribution of dark matter may 
be similarly complicated. 

3.2. Gravitational Lens Model 

Our goal is to constrain the shape of the mass distribution 
in the cluster core. We therefore adopt stringent criteria for 
the inclusion of multiple image systems as constraints on our 
lens model to guard against detection of spurious features in 
the mass distribution. To be included, a galaxy or morpholog- 
ical feature within a galaxy must be identified a minimum of 

among the California Institute of Technology, the University of California, 
and NASA. 



three times, and the morphological and color match between 
the multiple images of the galaxy/feature must be unambigu- 
ous. Ten morphological features of Al satisfy these criteria, 
of which eight are seen in all of Al.l, Al. 2 and A1.3, and two 
have also been identified in A1.4 (Fig. [TJ. These 10 features 
lie generally away from the portions of the disk that are af- 
fected by the cluster ellipticals marked with a black cross in 
Fig-Hi Note that we interpret all of the morphological features 
South-East of the central bulge (feature #2) of A 1.2 as being 
part of the disk of that image, distorted by the neighboring 
cluster elliptical. The three images each of A2, A3, A4, A6, 
and A7 are also used as model constraints, with the unknown 
redshifts of the latter three being free parameters in the lens 
model. These systems total n c = 61 model constraints. 

The mass distribution was initially parameterized as a su- 
perposition of 21 cluster galaxies (Ig\4 < 20.5), plus 5 cluster- 
scale components (hereafter referred to as halos) centered 
on the brightest galaxy in each of the light concentrations 
marked in Fig.Q] All galaxies and halos were parameterized 
as smoothly-truncated pseudo -isothermal elliptic al mass dis- 
tributions (PIEMD) following Ikneib et all {1996). The posi- 
tion, ellipticity, and orientation of the galaxies were matched 
to those of their light, and the velocity dispersions, core and 
cut-off radii were scaled with their luminosity, adopting the 
best-fi t parameters for an L* galaxy obtained by ISmith et all 
(120051) . We imposed a prior of 0.544 < z < 3 on the redshifts 
of A4, A6, and A7, the upper limit coming from the absence 
of an obvious Lyman break within/blueward of F555 W for all 
three galaxies. In total the model has n p = 16 free parameters. 

The model was fitted to the data usin g the Bayesian 
MCMC sampler within lenstool v6.5 12 dKneib et al.1 19961: 
Uullo et all2 007). using a positional uncertainty of 0.4" in the 
image plane. The image-plane x 2 of the best-fitting model 
was Xmin = 55.1, and the average root mean square (rms) 
deviation of images predicted by this model from the ob- 
served positions is (07) = 0.5". Halo A has a velocity dis- 
persion of a ~ 1270 kms -1 , and halos B, D, and E have 
a ~ 400 -500 kms -1 . However halo C has a 95% confidence 
upper limit of a < 343 km s" 1 , implying that the galaxies asso- 
ciated with halo C may not be embedded in an extended dark 
matter halo. 

To explore this further we re-fitted the model excluding halo 
C, again obtaining Xmin = 55-1 and (07) = 0.5". We then used 
the Bayesian evidence, i.e. the probability of the model given 
the data and the choice of the PIEMD paramaterization, to de- 
termine whether the additional complexity of the 5 halo model 
is justified by the data. The result is summarized in Table [T] 

- the probability of the 4 halo model exceeds that of the 5 
halo model by a factor of ~ 20 x . We therefore conclude that 
halo C is not justified by the data. We also test whether even 
simpler models offer more probable descriptions of the data 

- results are listed in Table Q] In summary, models with < 4 
halos are less probable than the four halo model by ~ 12—107 
orders of magnitude. We therefore adopt the four halo model 
as our fiducial model, and list its parameters in Table Q] 

3.3. Mass and Structure of the Cluster Core 

We use the fiducial model to compute the projected mass 
within a projected radius of R = 500 kpc and the fraction of 
that mass residing in the groups and the cluster galaxy pop- 
ulation, obtaining a mass of M(< 500 kpc) = (6.7 ±0.4) x 
10 14 M Q and a substructure fraction of / su b (< 500 kpc) = 

12 LENSTOOL is available online: http://www.oamp.fr/cosmology/lenstool 
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0.25 ±0.12. MACS Jl 149 therefore has a mass and substruc- 
ture fraction comparable with the m ost disturbed of the clus- 
ters stu died at z — 0.2 by Smith et al. (2005) and R ichard et alj 
(2009), see also Smith & Taylorl (120081) . However despite 
these global similarities, we note that none of the lower red- 
shift clusters contained three group-scale halos. MACS Jl 149 
is therefore the most complex strong-lensing cluster studied 
to date. 

Despite having similar velocity dispersions, the three 
group-scale halos have different optical morphologies, sug- 
gesting that they may have suffered different infall histories. 
Luminous structure B appears as an extended finger pointing 
SE from the cluster center (Fig. [TJ dominated by two bright 
elliptical galaxies: one on which we placed halo B in the lens 
model, and a second one ~ 20" further SE. Luminous struc- 
ture D comprises a dense group of galaxies, none of which 
dominate the optical luminosity. Finally, luminous structure 
E is dominated by a single bright elliptical galaxy. Despite our 
ability to find an acceptable fit, we therefore expect that im- 
provements to the strong-lensing constraints, especially spec- 
troscopic redshifts of A6 and A7 will help to constrain in more 
detail the structure of the se "groups", aided by more flexible 
modeling schemes (e.g. lJullo & Kne ib 2009). 

In Fig.|3]we show the contours of luminosity density Q3.1b , 
isomass density ( 33.21 Table [U, and X-ray surface b right- 
ness from Chandra observations (Ebeling et al. 2007). By 
construction the mass contours agree well with the luminos- 
ity density contours, except that luminous structure C is not 
embedded in an extended dark matter halo. The peak of the 
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FIG. 3.— The central ~ 1.5 Mpc xl.5 Mpc of MACS Jl 149, as revealed 
by the adaptively smoothed luminosity density of cluster galaxies (red), pro- 
jected total mass map calculated from the gravitational lens model (blue) and 
adaptively smoothed X-ray surface brightness contours from Chandra obser- 
vations (green). All contours are spaced linearly. 



X-ray emission is offset from the BCG by ~ 15", and the 
overall X-ray morphology is elongated in a NW-SE direc- 



Smith, G. P., et al. 



5 



tion, i.e. in the same direction as the mass and light con- 
tours. MACS J 1149 therefore follows the well-established 
trend for X-ray luminous clusters with multi-modal mass dis- 
tributions to have an X-ray morphology that is not centered 
on the BCG (ISmith et alJ2005HPoole et alJ2006trPowell et alJ 
2009; San derson et a l. 2009), adding weight to the conclusion 
that this is a merging cluster. 

3.4. Intrinsic Properties of the Spiral Galaxy at z=1.491 

We use the fiducial lens model (Table [U to calculate the 
intrinsic properties of the multiply-imaged disk-galaxy at z = 
1.491. Al.l, A1.2, AlJ, and A1.4 are magnified gravita- 
tionally by p, = 23, 18, 8, and 23 respectively, with a typical 
uncertainty of 30%. The unlensed apparent magnitude of this 
galaxy in the /-band is therefore / ~ 23 .4 ± 0.3, corresponding 
to Mb — -20.7 in the rest frame. Al is therefore comparabl e 
with L* galaxies in the local universe dNorberg et all [2 002), 
the fai ntest of the lensed disk galaxies in Swi nbank et alJ 's 
(2006) study of the Tully Fisher relati on at z — 1, and the 
lensed Sa galaxy at z = 1.6 studied bv lSmith etaT] d2002l) . In 
contrast, it is ~ 1 mag fainter than pr eviously studied unlensed 
disk galaxies at si milar redshifts (Ivan Dokkum & Stanford! 
1200 lH Wright et al.ll2008l) . Many of these galaxies are bulge 
dominated systems; in contrast, Al is not dominated by its 
bulge, with a rest frame B-band bulge-to-total ratio of B/T ~ 
0.4-0.5. 

We also estimate the star formation rate from the observe d 
V555-band flux - adopting the calibration of lKennicuttl (119981) . 
we obtain a global star formation rate of ~ 6 M Q yr -1 . Individ- 
ual HII regions in local galaxies ty pically span ~ 50- 100 pc 
dGonzalez Delgado & Perezll 19971) . which translates to ~ 6 - 
12 mas at z = 1.491. The lens magnification of p = 23 suf- 
fered by Al.l boosts the angular scale of HII regions in this 
galaxy to ^30-60 mas, bringing them within reach of instru- 
ments such as OSIRIS on Keck and NIFS on Gemini North. 
This galaxy therefore offers a unique opportunity to study the 
distribution of star formation, chemical abundance gradients, 
and even the physics of individual star-forming regions at a 
look back time of 9.3 Gyr at a level of detail similar to that 
achieved at z = 0.1. 

3.5. Density Profile 

We show the density profile from our fiducial model in 
Fig. S Parameterizing the profile as k ~ r 7 , we obtain 
7 ~ -0.3 ± 0.05 in t he radial range r ~ 3 — 30" . During peer 
review of this letter, Zitrin & Broa dhurstl d2009l hereafter ZB) 
claimed that the density profile of MACS J 1149 is flat and 
critically convergent (k ~ 1) out to r ~ 200 kpc, equivalent to 
~ 30". Inspection of ZB's Fig. 5 reveals that the average slope 
of their profile at r ~ 3 - 30" is 7 ~ -0.1 ± 0.02 (68% confi- 
dence, assuming their error bars are la). We therefore rule 
out ZB's model at ~ la confidence; a flat model (i.e. 7 = 0) 
within this radial range is ruled out at 12a. 

The two main differences between ZB's model and ours 
are the following. First, ZB's multiple-image interpretation 
is different from ours: they claim to find a fifth image of Al, 
plus six additional multiple-image systems (their 5-10), all of 
which do not pass the strict criteria described in 33.21 in addi- 
tion, ZB do not identify our systems A6 and A7. Second, ZB's 
analysis contains no spectroscopic or photometric redshift in- 
formation for their multiple images and, more importantly, 
ZB do not treat these unknown redshifts as free parameters 
in their model - their model contains just 6 free parameters, 
which describe the cluster mass distribution. 
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FIG. 4. — The density profile of MACS Jl 149 from our fiducial model 
(black solid), our ZB-constrained model with all redshifts as free parameters 
(red dashed), and the latter w ith th e redshifts of Al, A2, A3, A4 fixed (blue, 
dot-dashed), as described in j|3.5| In each case the grey filled regions show 
the 95% confidence interval around the best-fit model. The horizontal line 
marks the critical density required for strong lensing (k = 1 ), and the vertical 
dotted-dashed lines mark the average cluster-centric radius at which (from 
left to right) images of Al, A2, A3, and A4 are observed. 

We attempt to reproduce ZB's flat profile by fitting a model 
to all of their multiple image identifications, including the pu- 
tative fifth image of system Al (marked by a black circle in 
the A1.2 panel of Fig. 1). We treat the redshifts of all mul- 
tiple images as free parameters. The resulting best-fit "ZB- 
constrained" model has an image-plane rms of (er,) = 1.2" , 
i.e. more than twice that of our fiducial model, dominated by 
the fifth image of Al and ZB's systems 5-10. The density 
profile associated with this model is shown in Fig.|4]and is in 
fact steeper than ours. However, once redshifts are no longer 
included as free parameters in the fit, but fixed at values that 
differ, to varying degrees, from the true, measured values, the 
sensitivity of the density profile to a chosen set of fixed red- 
shift values becomes apparent. We demonstrate this by setting 
the redshifts of Al, A2, A3, and A4 to z = 1.5, 1.6, 1.8, and 
1.8, i.e. to values that are permitted by the model uncertain- 
ties but are in fact not the measured ones. The density profile 
(Fig. H]l resulting from these erroneous assumptions is nearly 
flat, with 7 ~ -0.14. 

We conclude that ZB's claim of a flat density profile 
is highly sensitive to the details of the method by which 
they chose to assign fixed redshifts to multiple-image sys- 
tems. These problems may have been compounded by mis- 
identification of some multiple image systems. 

4. SUMMARY 

We have presented new HST/ACS and Keck I/LRIS ob- 
servations of MACS J 1149, a massive X-ray selected galaxy 
cluster at z = 0.544 discovered in the Massive Cluster Survey. 
These data reveal seven robustly identified multiply-imaged 
galaxies, three of which we have confirmed spectroscopically. 
The most spectacular system is a multiply-imaged face-on 
disk galaxy at z = 1.491 that we identify as anL* (Mb ~-20.7) 
late-type (B/T <J 0.5) galaxy with an ongoing star formation 

rate of ~ 6 M yr" 1 ; the brightest images of this galaxy are 
magnified by p = 23. Future observations using integral field 
spectrographs should probe its properties in exquisite detail, 
thanks to the combination of lens magnification and fortuitous 
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viewing angle. 

We use the positions and redshifts of robustly identified 
multiply-imaged galaxies to constrain a detailed model of the 
mass and structure of the cluster core. Our fiducial model 
contains the main cluster halo plus three group-scale halos; 
the probability of a model this complex, relative to less com- 
plex models is P(iVhaio = 4)/P(/Y halo < 4) > 10 12 where AW> 
is the number of cluster/group-scale halos. We measure the 
mass and fraction of mass residing in substructures to be 
M(< 500 kpc) = 6.7 ± 0.4 x 10 14 M Q and / sub (< 500 kpc) = 
0.25 ±0.12 respectively. In summary, MACS Jl 149 is the 
most complex strong-lensing cluster core studied to date, its 
relatively dis-assembled nature being qualitatively consistent 
with the expectation that clusters at high redshifts are on av- 
erage less mature than those at lower redshifts. A more com- 
plete view will emerge from our analysis of the full sample of 
MACS clusters at z > 0.5 (Smith et al., in prep.). 

We also obtain a power law density profile slope of 7 = 
-0.3 ± 0.05 (95% confidence error bars) on scales of r ~ 3 - 
30", thereby ruling out density profile slopes as flat as those 
recently proposed by Zitri n & Broadhurst] d2009l) at J> la con- 



fidence. In summary. IZitrin & Broadhursl ' s result can be ex- 
plained by an absence of multiple-image redshifts of any form 
in their study, and by them not treating the unknown redshifts 
as free parameters in their model. These issues are probably 
compounded by them mis-identifying some multiple-image 
systems. Overall, this underlines the critical importance of 
measuring spectroscopic redshifts of multiply-imaged galax- 
ies for reliable lens models of strong lensing clusters. 
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